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ABSTRACT 

We revisit the occurrence condition of optically thick winds reported by iKato I (|1985f ) and 
iKato fc Hachisul (jT989) who examined mathematically nova envelope solutions with an old opacity 
and found that optically thick winds are accelerated only in massive white dwarfs (WDs) of > 0.9 Mq. 
With the OPAL opacity we find that the optically thick wind occurs for > 0.6 Mq WDs and that 
the occurrence of winds depends not only on the WD mass but also on the ignition mass. When the 
ignition mass is larger than a critical value, winds are suppressed by a density-inversion layer. Such a 
static solution can be realized in WDs of mass ~ 0.6 — 0.7 Mq. We propose that sequences consisting 
only of static solutions correspond to slow evolutions in symbiotic novae like PU Vul because PU Vul 
shows no indication of strong winds in a long-lasted flat peak followed by a very slow decline in its 
light curve. 

Subject headings: binaries: symbiotic — novae, cataclysmic variables — stars: individual (PU Vul) — 
stars: interiors — stars: mass loss 



1. INTRODUCTION 

Nova is a thermonuclear runaway event on a white 
dwarf (WD). During a nova outburst, the envelope on 
the WD expands to a giant size. It reaches the maxi- 
mum brightness and then becomes dark as the envelope 
loses its mass due to strong winds. Evolutions of such 
shell flashes have been theoretically followed by many 
authors with hydrodynamic/hydrostatic time-dependent 
calcu l ations (iPaczy hski & Zvtkow I [1 978: Sparks et al. 
19781: iNariaiet all a98Q: ilben i 1X982,: ,Prialnik et al. 



1981 iPrialnik fc Kovetz (ll99a iSion et al.lll979H bv us 

ing a sequence of steady-state solutions (jKato I [l983U 
IKato fc Hachisu ,1988. 1994), and by usmg a sequen ce 



of static solutions (jFuiimoto ifigsl iMacDonald Ill983f) . 

The time scale of nova evolution depends strongly on 
the WD mass. Massive WDs correspond to fast novae 
while less massive WDs do to slow novae. When the op- 
tically thick winds occur a large part of the envelope is 
blown off b y winds, and the nova d u ration is drastically 
shortened ()Kato fc Hachisu I 119941 : IPrialnik fc Kovetz I 
fl995l) . In less massive WDs winds are relatively weak, 
and in some cases, no optically thick winds occur. How- 
ever, the occurrence condition of winds has not been well 
st udied yet. 

iKato I (|1985f l mathematically examined the occurrence 
of optically thick winds with the Kramers opacity and 
found that the optically thick wind occurs when the 
luminosity approaches the Eddington luminosit y and 
the opacity increases outward. Kat o fc Hachisu I (|1989l ) 
showed that the optically thick wind occurs only in mas- 
sive WDs of MwD > 0.9 Mq, where Mwd is the WD 
mass. Wind acceleration is due to opacity enhance- 
ments at hydrogen and helium ionization zones. They 
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also showed that a nova evolution can be well repre- 
sented by a sequence consisting only of static solutions 
for Mwd < 0.9 Mq, and by a sequence of wind and static 
solutions for Mwd > 0.9 Mq. 

In the beginnin g of 1990's, opacity ta bles are revised 
(OPAL opacity: I glesias fc Rogers1ll996l and references 
therein), which show a prominent peak at logr(K)'^ 5.2 
strong enough to accelerate wi nds even in less massiv e 
WDs of Mwd ~ 0.5 - 0.6 Mq (IKato fc Hachisu I [1991 . 
However, the occurrence condition of winds has not been 
clarified yet, because the OPAL opacity has an intense 
peak at much higher temperature region than the helium 
ionization zone. Therefore, we re-examine the occurrence 
of optically thick winds for the OPAL opacity. We focus 
on less massive WDs in view of the application to slow 
novae or symbiotic novae in which less massive WDs cer- 
tainly inhabit. 

In [21 we introduce our simplified model for nova out- 
bursts. Then we present evolution sequences for three 
different WD masses in [JH Temporal changes of internal 
structure during nova outbursts are shown in [21 The 
occurrence condition of optically thick winds on various 
WD masses is presented in ij5l Discussion and summary 
follow in §Sj6]and|7l 

2. A SIMPLIFIED MODEL 

Figure [T] presents a schematic HR diagram of one cycle 
of nova outbursts. After thermonuclear runaway sets in 
on an accreting white dwarf (point A), the star bright- 
ens up and its envelope expands to a giant size. Optically 
thick winds blow in massive WDs during the extended 
envelope stage (from point B to point D through C). 
The star reaches the maximum radius at point C, where 
the envelope settles down into a thermal equilibrium in 
which nuclear energy generation is balanced with radia- 
tive loss. Afterwards, the star moves leftward keeping 
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Fig. 1. — A schematic HR diagram of one cycle of shell flashes. 
Thermonuclear runaway starts on an accreting WD at point A. The 
star brightens up and its envelope begins to expand. After the star 
reaches the maximum expansion at point C the star moves leftward 
along with the decreasing envelope mass sequence and the surface 
temperature increases with time. The nuclear burning extinguishes 
at point E and the star goes back to the original stage. When the 
optically thick wind occurs, it begins at point B and ends at point 
D. 

luminosity almost constant. Hydrogen nuclear burning 
stops at point E, and the star cools down to point A. 

When the wind occurs we solve the equations of mo- 
tion, mass continuity, radiative diffusion, and conserva- 
tion of energy, from the bottom of the hydrogen-rich en- 
velope through the photosphere assuming steady-state 
and spherical symmetry. This steady-state is a good ap- 
proximation in the decay phase or in w eak shell flashes 
i|Prialnik et al. 1 119861: [Kato et all 1 198 9^. When the op- 
tically thick wind does not occur, we solve the equation 
of hydrostatic balance instead of the equation of motion. 
In the rising phase, we integrated energy conservation 
equation without energy generation term due to nuclear 
burning and later estimated energy generation using the 
temperature and density obtained (Fujimoto 1982). In 
the decay phase we set the condition that the energy 
generation is balanced with radiative energy loss. These 
equation s and method of calcula tions are already pub- 
lished in ' Kato fc Hachisu I (|1994l ). 

In the rising phase, i.e., from point A to C in Figure 
[TJ we approximate the nova evolution by a sequence of 
envelope solutions as follows. When the optically thick 
wind occurs we use static solutions from point A to B 
and steady-state wind solutions from point B to C in 
Figure [H When the wind does not occur the entire se- 
quence consists only of static solutions. The occurrence 
of opticall y thick wind s is detected by the condition de- 
scribed in lKato I ()1985D :(1) the photospheric luminosity 
approaches the Eddington limit and (2) at the same time 
the thermal energy at the photosphere is comparable to 
the gravitational energy. We further assume that the en- 
velope mass does not decrease much in the rising phase, 
i.e., we approximate the rising-phase by a sequence of 
envelope solutions with a constant mass which we define 
as the "ignition mass", AAfig. This approximation may 
be too simple but we are interested in qualitative proper- 
ties of the envelope evolution, which are hardly changed 
even if we assume decreasing envelope mass. 



In the decay phase the envelope settles down into a 
thermal equilibrium, that is, the energy release due to 
nuclear burning is balanced to the outgoing energy flux. 
We approximate this stage by a sequence of optically 
thick wind solutions (from point C to point D) and static 
solutions (from point D to point E). When the wind does 
not occur, the entire decay phase are approximated by a 
sequence only of static solutions. In this sequence (from 
point C to E) the envelope mass is decreasing from its 
initial value of AAfjg at point C. Time-evolution is cal- 
culated from the mass decreasing rates by nuclear burn- 
ing a nd wind mass- loss if it occurs (see lKato fc Hachisu I 
[l99l for more details). Hydrogen burning stops at point 
E and after that, from point E to A, the evolution is fol- 
lowed by a static solution with a constant envelope mass. 

We have neglected effects of convection in wind so- 
lutions, because con vection is ineffective in a rapidly 
expanding envelope (|Kato fc Hachisu 1 11994[ ) . Convec- 
tive energy transport is calculated in static solutions 
using the mixing length theory with a = 1.5 other- 
wise specified. Effects of the a parameter on our re- 
sults are discussed in ^6.31 We use the OPAL opacity 
(jlglesias fc Rogers II1996I ). 

We further assume that the chemical composition of 
the envelope is uniform throughout the envelope and is 
unchanged with time, because the convection widely de- 
velops to mix the whole envelope in the early phase of 
outbursts. The solar composition is assumed for the en- 
velopes otherwise specified, which is a good approxima - 
tion for weak shell flashes. 'Prialni k fc Kovetzl (|1995l ) 
showed that the composition of cjecta is close to the so- 
lar for high mass accretion rates. 

The above approximations may not be accurate enough 
to follow nova outbursts but sufficient for our purpose of 
qualitative study on the occurrence of winds and on the 
internal structures of WD envelopes. 

3. NOVA SEQUENCES 

Figures [H [3l and |4] show the quasi-evolution sequences 
that mimic the rising and decay phases of nova out- 
bursts. In a case of 0.5 Mq WD, no optically thick 
winds are accelerated at all, thus all the sequences in 
Figures [2] consist of static solutions. Figure [2^ shows a 
case with a small ignition mass of AMig = 4.0x10"^ Mq . 
The rising phase ends when the photospheric tempera- 
ture and radius reach logT (K)= 5.25 and logr (cm)= 
9.77 (0.085 Rq), respectively. This point is the maximum 
expansion where the envelope reaches thermal equilib- 
rium. The maximum expansion shifts toward lower tem- 
perature when we increase the ignition mass. In a case of 
AAfig = 5.0 X lO"'^ Mq (Fig. Eb) the temperature and 
radius at the maximum expansion are logT (K)= 4.8 
and logr (cm)= 10.70 (0.73 Rq), respectively, and for 
AMig = 1. X 10--* Mq (Fig. Eb), logT (K)= 3.79 and 
logr (cni)= 12.71 (74 Rq), respectively. 

Figure [3] shows a 1.0 Mq WD case. If we increase 
the ignition mass as AMjg = 3 x 10"'', 8 x 10"'^ Mq, 
and 4 x 10~^ Mq, the temperature and radius at the 
maximum expansion change to logT (K)= 5.6, 4.7, and 
3.8, and logr (cm)= 9.4 (O.O36i?0), 11.18 {2.2Rq), and 
12.83 (97i?0), respectively. We found that optically thick 
winds occur for AA/jg > 3.6 x 10"^ Mq. Wc define 
this critical mass for winds as AAf^ind, i-e., winds oc- 
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Fig. 2. — Rising (black) and decay (red) phases of a nova outburst 
on a 0.5 Mq WD. The position of maximum expansion is indicated 
by a small open circle, (a) AAfig = 4 X 10~^ Mq, (b) 5 X 10~^ Mq, 
and (c) 1 X 10"'' Mq . 
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Fig. 3.— Same as Fig[2l but for a 1.0 Mq WD. (a) AMig = 
3 X lO"** Mq, (b) 8 X lO"*^ Mq, and (c) 4 X lO'^ Mq. The solid 
and dashed lines indicate static and wind phases, respectively. The 
position where the wind occurs/stops is denoted by a cross. 
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Fig. 4.— Same as FigO but for a 0.6 Mq WD. (a) AMig = 
2.2 X 10-5 Mq, (b) 4.2 X 10"^ Mq, and (c) 7 x 10"^ Mq. 

due to the opacity peak at logT (K)'^ 5.2, optically 
thick winds occur in the lowe r temperature side of the 
peak (|Kato fc Hachisu I Il994f ) as shown by the dashed 
line in Figures and [3J;. If we further increase the 
ignition mass to AMig > 8.5 x 10^"* Mq, winds are sup- 
pressed and static solutions exist instead of wind solu- 
tions. This is because hydrostatic balance is estabhshed 
in the envelope when a substantial amount of matter 
above the super- Eddington region, i.e., around the opac- 
ity peak, supresses down the wind-acceleration. We de- 
fine this critical mass for static expansion, as AMoxp- If 
AAfig > AMoxp no winds are accelerated and the en- 
velope expands without optically thick winds. However, 
this critical value is too l arge for the ignition mas s of 
an accreting 1.0 Mq WD (jPrialnik fc Kovetz I [19951 ). so 
such static solutions are not realized and we don't go into 
details. 

These two critical masses of AMwind and AMexp de- 
pend on the WD mass. For a 0.6 Mq WD, we obtained 



AM„ 



2.3 X 10"^ Mq and AM, 



6.5 X 10-5 Mq 



Figured^ shows a case of AMjg = 2.2 x 10"^ Mq, in 
which no winds occur because the photospheric temper- 
ature does not reach the OPAL peak even at the max- 
imum expansion at logT (K)= 5.35 and logr (cm)= 
9.68 (O.O69i?0). In the case of AMig = 4.2 x 10^^ Mq 
(Fig. [JId) the optically thick winds begin when the 
photospheric temperature decreases to logr(K)= 5.25. 
The maximum expansion reaches the photospheric ra- 
dius of logr(cm)= 11.60 {5.7Rq) and the temperature of 
logT(K) = 4.39. Figure it shows the case of AMig = 
7.0 X 10-5 Mq (> AMoxp) in which the whole period 
of the outburst can be followed by a sequence of static 
solutions. 
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Fig. 5. — Evolutional change of the diffusive and Eddington 
luminosities in the rising phase of Figure |4j) (the envelope mass of 
AAfig = 4.2 X 10-5 Mq on the 0.6 Mq WD). The top three panels 
show the static solutions for different I, the ratio of the luminosity 
to the Eddington luminosity at the photosphere. The bottom three 
panels show the optically thick wind solutions in which the small 
dot denotes the critical point l[Kato 1983a). The right edge of 
each curve corresponds to the photosphere, the radius of which is 
(a) logi?ph (cm)=9.08, (b) 9.23, (c) 9.86, (d) 10.05, (e) 11.33, and 
(f) 11.60. The model (a) locates under the "knee" of the line in 
Fig. |4]3 so the photospheric radius and temperature are smaller 
than those of model (b). 

Note that there is only one path for the decay phase of 
the 0.5 Mq WD when it is in a thermal equilibrium. The 
two decay phases in Figures [2^ andEt" are a part of that 
in Figure Similarly the decay phases of Figures [3^ 
and [Sb are identical to that of Figure [S};. On the other 
hand, there are two different paths in the 0.6 Mq WD 
even when it is in thermal equilibrium. The decay phase 
of Figure |4)3 consists of wind and static solutions, only 
the latter part of which is identical to that of Figure Hj;. 
The wind sequence in Figure |4|d and the static sequence 
at logTph < 5.29 in Figure [4); are different from each 
other all in envelope structure, in evolution time scale, 
and in light curve as we will see later. 

4. INTERNAL STRUCTURES 

Figure [5] shows the distribution of the diffusive lumi- 
nosity and the local Eddington luminosity against the 
temperature for solutions along the rising phase of Fig- 
ureHb, i.e., 0.6 M© WD of AA/jg = 4.2 xlO'^ Mq. Here, 
the local Eddington luminosity is defined as 



^Edd = 



AttcGM 



(1) 




log 



(cm) 



where k is the opacity in which we use the OPAL opacity 
(jlglesias fc Rogers Ill996l ). Since the opacity k is a func- 



FlG. 6. — Evolutional change of internal structure of the en- 
velope in the rising phase of Figure (0.6 Mq WD with 
AMig = 4.2 X 10~^ Mq). From top to bottom, the wind velocity 
(V), temperature (T), density (p), Eddington luminosity (Z/Edd) 
and diffusive luminosity (Lt). The dashed line denotes the escape 
velocity defined by 2G MYfn / r . The right edge of each curve 
corresponds to the photosphere. Dot denotes the critical point of 
each wind solution. 



tion of the temperature and density, the Eddington lumi- 
nosity is also a local variable. This Eddington luminosity 
has a local minimum at logT (K) = 5.25 corresponding 
to the opacity peak as shown in Figure [Bf. 

As the star moves upward in Figure |4)d, the diffusive 
luminosity increases and approaches the Eddington lu- 
minosity near the photosphere (Figures [n^,b, and c). 
When the photospheric temperature decreases to logT 
(K) ^ 5.2, matter is accelerated and steady mass-loss be- 
gins (Figure [5ji). After that, the envelope continuously 
expands to until point C in Figure [TJ A narrow super- 
Eddington region appears corresponding to the opacity 
peak at logT (K) ~ 5.2. 

The change of envelope structure is shown in Figure 
m We see that the structure of the static solution just 
before the wind occurs is very similar to that of the ad- 
jacent wind s olution . This property has been already 
pointed out bv lKato I ((1985. ) for the old opacity and here 
we confirm it for the OPAL opacity having a prominent 
peak at a much higher temperature region. 

In the massive envelope of AMig = 7.0 x 10^^ Mq, 
however, no winds arise when the envelope expands be- 
yond the opacity peak of logT (K) ~ 5.2. Figures [7| and 
[8] show internal structures in such an evolution sequence. 
Difference from the wind sequence becomes prominent as 
the envelope expands (see Figs. [5] and ^ . 

A remarkable difference can be observed in the den- 
sity distribution. Figure [5] shows the monotonically de- 
creasing density having a dependence in the outer 
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Fig. 7. — Same as Fig. [5] but for the rising phase of Fig. |4j; 
(0.6 Mq WD with AA/ig = 7.0 x 10-^ Mq). The photospheric 
radius of each solution is (a) logij (cm)=9.11, (b) 9.50, (c) 9.96, 
(d) 10.56, (e) 11.16, and (f) 12.37. 

envelope (logr (cm) > 10.3) as expected from the equa- 
tion of continuity Anr"^ pv = const, where the velocity 
V is almost constant in an outer part of the envelope. 
On the other hand, the envelope in the static sequence 
(Fig. [8|) develops large density-inversion at logr (cm) 
~ 10 — 11.3 corresponding to a super Eddington region. 
This density-inversion arises in order to keep hydrostatic 
balance in the super-Eddington region (LEdd < Lr) as 
expected from the equation of hydrostatic balance. In- 
efficient convections occur in the region of ^Edd < Lr 
but are unable to carry all of the diffusive energy flux, 
thus the structure is super-adiabatic. Comparing Figure 
[8] with Figure [6] we see a prominent core-halo structure 
in the density and temperature distributions in the static 
solutions. 

Such differences between the wind and static sequences 
can be also seen in the decay phase. Figure [S] demon- 
strates the difference between the two solutions in the 
decay phase of 0.6 Mq WD. Both the static and wind so- 
lutions have the same photospheric temperature logTph 
(K) ~ 4.53 but they are very different in their internal 
structures. The mass-losing envelope shows the density 
distribution of r^^ in the outer part while the quasi-static 
envelope develops a wide density-inversion region. 

5. OCCURRENCE OF OPTICALLY THICK WINDS 

As we have seen in the previous sections, the wind 
mass-loss occurs in a limited range of the envelope mass. 
When the ignition mass is smaller than AM^ind the 
envelope expands a little and the photospheric tem- 
perature does not reach the opacity peak of logT (K) 
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Fig. 8. — Same as Fig. |6] but for the rising phase of Fig. \^ (0.6 



Mq WD with AMig = 7.0 x 10"^ Mq). 
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Comparison of two solutions with the same photo- 



spheric temperature of logTpjj = 4.53 in the decay phase of 0.6 
A/0 WD. Thick line: envelope solution in the static sequence. The 
envelope mass is 2.3 X 10~^ Mq. Thin line: envelope solution in 
the wind sequence. The envelope mass is 3.8 X 10~^ Mq. The 
right edge of each line corresponds to the photosphere. 



^ 5.2 (see Figs. ^jji. and 2^). Therefore, no optically 
thick winds occur. On the other hand, if the ignition 
mass is larger than AMexp, winds are suppressed in a 
way that density-inversion balances to radiation-pressure 
gradient in a super-adiabatic region (see Fig. [8]). In 
this case the envelope expands without optically thick 
winds. Therefore, optically thick winds occur only for 
AMwind < AMig < AA/cxp as we have already seen in 
the 0.6 Mq WD. 

Figure [TU] depicts these two critical masses of AAf„ind 
and AMoxp for various WD masses. Optically thick 
winds are driven in the right hand side of the thick solid 
line. In the lower region (AMjg < AM^ind) the enve- 
lope does not expand much and no optically thick winds 
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Fig. 10. — Upper and lower critical ignition masses for winds, 
AMcxp and AAf^inj, plotted against the WD mass. The opti- 
cally thick winds occur in the right hand side of the lines (la- 
beled "wind"). In the upper side of the line, optically thick wind 
does not occur and the envelope expands quasi-statically (labeled 
"expansion"). In the lower region, the envelope does not ex- 
pand so much ("no expansion") and no wind arises. Thick line: 
solar composition. Dashed line: CNO rich {X = 0.35, Y = 
0.33, XcNO = 0.30,and Z = 0.02). Dash-Dotted line: Popula- 
tion II composition {X = 0.7, Z = 0.004). Open circles with a 
dot an d crosses denote mass ejec tion model and no mass ejection 
model llPrialnik fc Kovetz 1119951) . respectively. See text for more 
details. 

arise. In the upper region (AMjg > AMoxp) the envelope 
is so massive and winds are suppressed. Here, we label 
the lower region "no expansion" , and the upper region 
"expansion" , which is an abbreviation of "quasi-static 
expansion" . 

As the OPAL opacity depends on the chemical compo- 
sition of the envelope, both AMgxp and AMwind depend 
on the composition. For a composition of typical classical 
novae {X = 0.35, F = 0.33, Xcno = 0.30, and Z = 0.02) 
the wind is strongly accelerated and the "wind" region 
extends as shown in Figure [ini On the other hand, for 
Popul ation II stars with lower metal content, the wind is 
weak (jKato lll997lll999D and the "wind" region becomes 
narrower. 

Figure [10] also shows that AMexp increases with the 
WD mass and reaches as large as 10"'^ — 10~^ Mq for 
1.0—1.38 Mq WDs. Such a large ignition mass is unlikely 
to be realized in accreting WDs. Dynamical calculations 
have shown that the envelope mass is up to 2 x lO^'' Mq 
for a 1.0 Mq WD and 3 x 10"'^ ^.^^ fo^. ^ ^ 4 foj. ^^e 
accretion r ate of 10"^^ '^ Mq y r~^ and 10~^^ Mq yr~^ for 
cold WDs (lYaron et al.l 120051 ). These values are practi- 
cally upper limits for the envelope mass in nova outbursts 
but are still much smaller than AMexp- Therefore, we re- 
gard that the " expansion" region and a part of the upper 
"wind" region theoretically exist but may not be realized 
in the actual mass-accreting 1.0 — 1.38 Mq WDs. 

In less massive WDs such as ~ 0.6 Mq, on the other 
hand, AAfcxp is as small as 10~^ Mq, which corre- 
sponds to the mass-accretion rates of ^ 10^^ Mq yr~^ 




t (day) 

Fig. 11. — Comparison of light curves in the decay phase 
of a 0.6 Mq WD between the "wind" sequences (dashed) and 
"static" expansion sequences (solid) for the same chemical compo- 
sition. Dotted line: Static expansion sequence with mixing length- 
parameter a = 2.0. Solid line: a = 1.5. Dash-dotted line: a = 1.2. 

(jPrialnik fc Kovetz II1995D . a typical mass accretion rate 
in cataclysmic variables. Therefore, the "expansion" 
region becomes a subject of realistic interest. If the 
ignition mass is larger than Ailfoxp the envelope ex- 
pands in a quasi-static manner. In less massive WDs 
of AfwD < 0.5 Mq the "wind" region disappears com- 
pletely and the envelope evolves in a quasi-static manner. 

6. DISCUSSION 
6.1. Light Curves in the Decay Phase 

As we see in Figure |4] there are two different sequences 
that represent the decay phase of nova outbursts on a 
0.6 Mq WD. These two sequences are different from each 
other in its envelope mass and mass decreasing rate, so 
the light curves are also different. Figure fTD shows the 
light curves corresponding to these sequences, both of 
which start from logTph (K) ~ 3.9. The light curve of 
the static sequence decays much slowly, because its evolu- 
tion speed is determined by the mass-decreasing rate due 
only to hydrogen nuclear burning, whereas in the "wind" 
sequence the evolution is accelerated by the winds which 
carry away a large part of the envelope matter. The most 
remarkable difference is in the flat peak of the static se- 
quence, whereas it decays sharply in the wind sequence. 

The flat peak in the static sequence reminds us 
a peculiar light curve of the symbiotic nova PU 
Vul, which shows a very long plateau peak of 3000 
days followed by a slow decline of 3 magnitudes/2500 
days. The spectra suggest a very quiet expansion 
with no indication of wind rnass loss during the fiat 
peak (lliiima fc Ortolani II1984I: iKanamitsu et al. 1 119911 : 
iKolotilov et al.lll995l: lYamashita et al. II1982D . We may 
apply our static "expansion" sequences to such slow out- 
bursts. Detailed light curve fitting will be presented sep- 
arately. 

Both of the wind and static expansion solutions are cer- 
tainly stable because they represent nova envelope and 
red giant envelope, respectively. However, it is interest- 
ing to point out a possibility that a nova evolves along 
the static expansion sequence but it suddenly jumps in 
the wind sequence or vice versa during the course of 
evolution. If this kind of transition occurs it may pro- 
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ceed from a higher total enthalpy state to a lower one. 
The enthalpy integrated for the entire envelope is, for 
example, 3.6 x 10^^ erg for a static envelope solution 
of AM = 5.4 X lO-'^ Mq with logTph (K) = 3.8, and 
1.0 X 10^^ erg for a wind mass-loss solution of a similar 
envelope mass and photospheric temperature. Therefore, 
if the transition occurs it possibly goes from the static 
expansion sequence to the wind sequence. As the inter- 
nal structure of the two solutions are very different (see 
Fig. [HI) and the excess energy will be released, the tran- 
sition may not occur in a quiet way but accompany some 
violent activities. It may be interesting to connect such 
activities to oscillatory behaviors often observed in early 
light curves of slow/symbiotic novae. 

Figure [Til demonstrates that both of the static "expan- 
sion" and "wind" sequences exist for the same WD mass 
and the same chemical composition of the envelope. If 
the accretion rate onto the WD changes with time and, 
as a result, the ignition mass changes from one outburst 
to the next around AA/cxp, these two different types of 
outbursts can be realized on the same WD. In such a 
case outbursts behave very differently for a small change 
of mass-accretion rate. 

In this way, we expect many active phenomena associ- 
ated with static expansion sequences. More quantitative 
studies including light curve analysis of slow/symbiotic 
novae are necessary. 

6.2. Comparison with Dynamical Calculations 

We have obtained AM^ind, the lower critical mass for 
the envelope having optically thick winds. This criti- 
cal inass_2S_coni£ajedwitl^ calculations 
by [Prialnik fc Kovetzl (|1995f ) who presented a number 
of shell flashes for various parameters. The open circles 
with a dot in Figure [TOl denote their models in which 
mass ejection occurs, and the crosses indicate the mod- 
els with no mass ejection. In these models the chemi- 
cal composition of the envelope is not uniform but close 
to the solar value because dredge-up of WD material is 
insignificant for high mass-accretion rates. Therefore, 
these models can be compared with our model labeled 
"solar" . Our critical line is quite consistent with their 
results except the 0.65 Mq WD, in which their crosses 
are above our line. These authors defined the "mass ejec- 
tion" by that the expansion velocity exceeds the escape 
velocity, which is different from our definition of occur- 
rence condition of the winds described in 2) In our wind 
solutions of 0.6 M© the velocity does not exceeds the es- 
cape velocity as shown in Figure [6] Considering this 
difference and other difference of input parameters, we 
ca n say that our AMymd is quite consistent with results 
in [Prialnik fc Kovetzl ([1995[ ). 

6.3. Dependence on the Mixing Length Parameter 
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We assumed a mixing-length parameter of a = 1.5. 
Many authors have estimat e d the a pararneter t o be 
a = 1.2 - 2.0 ([Asida I [2OOOI : [Palmieri et aTl [20021 and 
references therein) for various types of stars. The mixing- 
length parameter could be a function of stellar param- 
eters or could depend on stellar structure, but these 
dependences are not known yet. Therefore, we simply 
adopted a — 1.5. 

In order to see the dependence of light curves on the 
mixing-length parameter we have calculated additional 
models with a = 1.2 and 2.0 as shown in Figure [TTl For 
a smaller a, the energy transport is more efficient and 
the star evolves faster, resulting in a steeper light curve. 
Even though, the static sequences still show a long-lasted 
flat peak and slow decline after that. 

7. SUMMARY 

Our main results are summarized as follows; 

(1) For a given WD mass and chemical composition, 
the optically thick wind occurs when the ignition mass 
(AMig) satisfies the condition, i.e., AM^ind < AMig < 
AMcxp- If AMig < AMwind the envelope does not ex- 
pand and no wind is accelerated. When AMig > AMexp, 
winds are suppressed and the envelope expands with no 
optically thick winds. 

(2) Optically thick winds occur smoothly from a static 
envelope because the structure is not drastically changed 
before and after t he win d occu rs. This property was 
already reported in [Kato I ([1985[ ) for the old opacity but 
we confirm this for the OPAL opacity. 

(3) For a given WD mass and chemical composition, 
there exist two solutions with different structures. One 
is the wind solution with a monotonic decrease in the 
density distribution. Another is the static solution that 
develops a wide super-adiabatic region with density in- 
version. 

(4) In massive WDs (> 0.8 Mq), AMcxp is very large 
and the static-expansion sequence is not practically real- 
ized (AMig «C AMexp). In the intermediate mass WDs 
(0.6 — 0.7 Mq), both the optically thick wind and quasi- 
static expansion are realized depending on the ignition 
mass. In less massive WDs (< 0.5 M©) no optically thick 
wind occurs independently of the ignition mass. 

(5) The wind sequence has been applied to nova out- 
bursts in which strong mass-loss is observed. The newly 
found quasi-static expansion sequence may be applied 
to slow evolutions of symbiotic novae like PU Vul, which 
shows a long-lasted flat peak with no indication of strong 
winds. 
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